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The mono-jet search, looking for events involving missing transverse energy (j5 T ) plus one or two 
jets, is the most prominent collider dark matter search. We show that multi-jet searches, which 
look for $ T plus two or more jets, are significantly more sensitive than the mono-jet search for 
pseudoscalar- and scalar-mediated interactions. We demonstrate this in the context of a simplified 
model with a pseudoscalar interaction that explains the excess in GeV energy gamma rays observed 
by the Fermi Large Area Telescope. We show that multi-jet searches already constrain a pseudoscalar 
interpretation of the excess in much of the parameter space where the mass of the mediator Ma 
is more than twice the dark matter mass judm. With the forthcoming run of the LHC at higher 
energies, the remaining regions of the parameter space where Ma > 2ttidm will be fully explored. 
Furthermore, we highlight the importance of complementing the mono-jet final state with multi-jet 
final states to maximise the sensitivity of the search for the production of dark matter at colliders. 


INTRODUCTION 

Weakly interacting massive particles (WIMPs) are the 
most studied and arguably the best motivated candidate 
for particle dark matter (DM) as they are present in many 
extensions of the Standard Model (SM). A particularly 
appealing feature of WIMPs is that they should be de¬ 
tectable with current or near-term experiments [f]. 

The plethora of DM models poses a challenge of how to 
interpret DM searches in a generic way. One approach 
is to classify the DM model by the particle mediating 
the interaction. A particularly interesting class of mod¬ 
els involves the exchange of a spin-0 s-channel scalar or 
pseudoscalar mediator, since additional scalars and pseu¬ 
doscalars are a generic prediction of extensions of the 
Standard Model (SM) Higgs sector. Pseudoscalars are 
also particularly interesting as they are a common feature 
in many of the models proposed to explain the spatially 
extended gamnra-ray excess around the Galactic Centre 
observed with the Fermi Large Area Telescope (Fermi- 
LAT) [HE]. 

While scalars and pseudoscalars with a mass below 
10 GeV can be probed by flavour-changing observables 
at colliders |4], heavier pseudoscalars whose dominant 
interaction is with DM are particularly difficult to de¬ 
tect. Pseudoscalar-mediated interactions result in a sup¬ 
pressed tree-level spin-dependent interaction and an un- 
observably small loop-level spin-independent interaction 
at direct detection experiments, making this interaction 
inaccessible for these experiments W- The most promi¬ 
nent collider search for DM production at the Large 
Hadron Collider (LHC) is the mono-jet search [TUTU], 
which searches for events with a high momentum jet from 
initial state radiation in combination with significant 
missing transverse energy ($t)- As we will demonstrate 
(see also fTTUISl ). the mono-jet search has limited sensi¬ 
tivity to pseudoscalar- and scalar-mediated interactions. 


In contrast, we show that multi-jet plus lfc T collider 
searches significantly extend the sensitivity of the LHC 
to these interactions. These searches are designed to 
be inclusive and probe a large region of the topological 
and kinematic phase space, probing jet-multiplicities > 2 
with several kinematic variables, including ljl T and the 
scalar sum of the jets pp ( Ht )• Typically, the multi-jet 
plus $ T final state has been used to search for supersym¬ 
metry (SUSY) at the LHC. In this letter we demonstrate 
that this final state also has excellent sensitivity to the 
pair-production of DM from pseudoscalar and scalar me¬ 
diators. This is because the production of pseudoscalar- 
or scalar-mediators is typically dominated by gluon fu¬ 
sion m, which in turn generally leads to events with 
higher jet-multiplicity in the final state [T5j. 

As an example of the utility of the multi-jet plus 1$ T 
searches, we apply our limits to a pseudoscalar model 
that can account for the Fermi-LAT excess. We show 
that current limits already exclude much of the Fermi- 
LAT excess parameter space where the pseudoscalar’s 
mass (Ma) is more than twice the DM mass (todm)- Fu¬ 
ture limits at y/s =13 TeV will fully probe the remaining 
parameter space where Ma > 2wdm- 

MONO-JET AND MULTI-JET SEARCHES 

The benchmark pseudoscalar model that we consider 
is a simplified model following the ansatz of the Minimal 
Simplified Dark Matter (MSDM) models T6], which have 
four free parameters: tttdm, Ma and two couplings guM 
and gsM- The interaction terms are 

Ant = ignMA-xfx + *3sm V — A q~f>q , (1) 

Z —' V 

Q 

where x is a Dirac fermion, the sum is over all quarks, m q 
is the quark mass and v = 246 GeV is the Higgs vacuum 
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expectation value. Motivated by the Minimum Flavour 
Violation hypothesis HZ3, we assume that the couplings 
of the pseudoscalar to quarks are proportional to m q . 
With this coupling structure, the s-channel production 
of A is dominated by gluon fusion. We adopt a simpli¬ 
fied model approach in this work as it provides a more 
accurate framework to characterise the results of collider 
DM searches when the mediator is light enough to be 
produced mmm- This simple model can explain the 
Fermi-LAT excess while remaining consistent with other 
constraints jBJ, and is a useful proxy for the structure 
found in two-Higgs-doublet models (2HDM) and in ex¬ 
tended 2HDM that have mixing with a singlet-like pseu¬ 
doscalar (a la the NMSSM) [T9] . 

The MSDM ansatz assumes that the pseudoscalar 
width T a is fully determined by the decays to quarks, 
gluons and DM. The only free parameters affecting 
are the four basic parameters {todm, M a , 5dm, 5sm}- Ex¬ 
pressions for T a are given in Ref. m- r A is domi¬ 
nated by DM and top-quark decay (when kinematically 
allowed) as m q /v does not suppress these decays. We 
will focus on the regime where M A > 2?txdm in which 
case BR(A -A \x) is generally large [20]. The collider 
signatures therefore involve the production of A that de¬ 
cays to a pair of DM particles. 

The mono-jet search at the LHC has been advocated 
as the primary model-independent search for DM pro¬ 
duction at colliders [3 IS]. Events with a second jet may 
also be allowed but events with three or more are rejected 
to avoid background contamination from processes with 
high jet-multiplicity, like top-quark production. 

In contrast, the more inclusive cut [H|, Mt 2 [22], Ra¬ 
zor [23], or MHT-Ht [21] multi-jet plus $ T searches 
place fewer constraints on the phase space. Each event 
is characterised by the number of jets and hadronic ac¬ 
tivity Ht, as well as other kinematic variables. These 
bins are combined in a likelihood fit and allow the multi¬ 
jet searches to take advantage of different signal-to- 
background compositions in these numerous search re¬ 
gions to attain better sensitivity. For instance, whilst 
the CMS mono-jet analysis employs a single inclusive $ T 
bin with the threshold ranging between 250 GeV and 
550 GeV, the M ?2 search combines more than 100 ex¬ 
clusive search regions. Similar kinematic selections and 
jet-multiplicity categorisations are utilised by the cut and 
Razor searches. These inclusive searches are an impor¬ 
tant pillar of the search strategy for new physics at the 
LHC, providing the best possible sensitivity to a large 
variety of SUSY production and decay topologies [25] , 
So far they have largely been ignored for searches in¬ 
volving the pair production of DM. A Razor search was 
previously investigated in [26] [27] but found a small im¬ 
provement over the mono-jet search. The search in [26] 
considered one inclusive signal region, which was opti¬ 
mised for models with a (axial-)vector mediator. Most of 
the events arising from gluon-induced models, like scalar 


or pseudoscalar exchanges, were rejected. In contrast, 
our analysis includes all of the accessible signal regions 
of the Mt2 search and maintains excellent performance 
for a variety of signal models. Unfortunately, Ref. m 
did not consider gluon-induced models so no direct com¬ 
parison is possible. 

To determine the sensitivity of these collider searches 
for our model, we reinterpret the CMS mono-jet [2] and 
Mt2 [221 analyses using the Powheg Box V2 genera¬ 
tor [15] [28] . This generates, at leading order with ex¬ 
act top-quark mass effects, DM pair production together 
with one parton via an s-channel pseudoscalar mediator. 
We use the fixed width approximation, having checked 
that our results match when the running width is used. 
We use the MSTw2008lo parton distribution functions 
with renormalisation (/x^) and factorisation scale (/i_p) 

set to 5/2, where /x = + P? ji +Ptj 1, m xx the 

invariant mass of the DM pair and PT,ji is the transverse 
momentum of the leading jet. Scale uncertainties on the 
cross-section were found to be 0(±4O%) [13] . To be con¬ 
servative, we do not apply a A'-factor of 1.6 as used in EH 
to account for higher order corrections, since a computa¬ 
tion of the next-to-leading-order corrections with a finite 
top-quark mass is not available. The events generated 
by Powheg are interfaced with Pythia 8.180 [22] for 
parton-shower effects and hadronisation. Finally, signal 
events are passed through Delphes V3.2.0 [30] for de¬ 
tector simulation. 

Figure [l] shows a comparison of the expected 90% con¬ 
fidence level (CL) exclusion contours from our mono¬ 
jet and Mt2 analyses in the M A — mo m plane (upper) 
for two different coupling scenarios, 5sm = 5dm = 1 
and 3, as well as in the 5sm — 5dm plane (lower) for 
wdm = 45 GeV and M A = 150 GeV and 250 GeV. 
This comparison assumes 20 fb -1 at 8 TeV. Both planes 
show that the Mt 2 search possesses significantly better 
expected sensitivity over the mono-jet search. In par¬ 
ticular, the Mt2 search has the ability to exclude large 
regions in the M A — ?tt,dm plane for both 5sm = 5dm = 1 
and 3, up to M A = 350 GeV and 550 GeV respectively 
for small ttidm- In contrast, for 5sm = 5dm = 1 the 
mono-jet analysis does not find any limit while the reach 
in M a for 5sm = 5dm = 3 is more than 150 GeV less than 
for Mt 2- The lower panel shows that for M A > 2?bdm, 
the multi-jet plus $ T analysis is expected to probe cou¬ 
plings down to 5dm — 0.05 and 5sm — 1- In comparison, 
the mono-jet search possesses sensitivity for 5 dm > 0.12 
and 5 sm 2, consistent with the findings in [13] . 

Based on these results, we conclude that multi-jet 
plus lfi T searches exhibit better sensitivity than the 
mono-jet analysis over the entire parameter space. The 
improved sensitivity of Mt2 is a result of categorising the 
search in bins of jet-multiplicity and several kinematic 
variables, allowing for differences in signal and back¬ 
ground in these various categories to be exploited. We 
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FIG. 1 . Comparison of the expected 90 % CL exclusion con¬ 
tours from our mono-jet (blue line) and Mt2 (red line) anal¬ 
yses. Regions below and above the lines are excluded in the 
upper and lower panels, respectively. The expected Mt2 lim¬ 
its are significantly better than the expected mono-jet limits 
over the entire parameter space. 

find that the low and medium Ht categories for two jets 
provide a large fraction of the sensitivity for our pseu¬ 
doscalar model. However, significant additional sensitiv¬ 
ity is gained by the inclusion of low and medium Ht cat¬ 
egories with 3-5 jets. The higher-jet bins are particularly 
important for our model since ~ 60% and ~ 30% of the 
events in the low and medium Ht categories have gluon- 
fusion (gg) and quark-gluon (qg) production, which typi¬ 
cally produce more jets in the final state. The remaining 
~ 10% of events are from quark anti-quark (qq) or gluon 
heavy-quark initial states. Having demonstrated the en¬ 
hanced sensitivity of Mt 2 over the mono-jet search, in the 
following we will show only limits from the Mt 2 analysis. 


CONSTRAINING THE FERMI-LAT EXCESS 

The Fermi-LAT observation of a spatially extended 
gamma-ray excess around the Galactic Centre has gener¬ 


ated great interest since it may be explained by annihilat¬ 
ing DM. Unfortunately, many indirect detection signals, 
including the Fermi-LAT excess, do not give irrefutable 
evidence for DM because of large astrophysical uncertain¬ 
ties m- For instance, Ref. J32] suggests that the excess 
could be explained by point sources (PS) that lie just be¬ 
low the current Fermi-LAT threshold. While detecting 
members of the PS population would corroborate an as¬ 
trophysical origin for the excess, a complementary signal 
in direct detection or collider experiments is required to 
corroborate a DM origin. 

A plethora of models involving a pseudoscalar me¬ 
diator have been proposed to explain the Fermi-LAT 
excess O [33]. As pseudoscalar-mediated interactions 
are suppressed at direct detection experiments, collid¬ 
ers are the most promising way to independently test 
a pseudoscalar-mediated explanation for the gamma-ray 
excess. We therefore investigate the implications of the 
Mt 2 limits on the model defined by Eq. ([l]), which can 
explain the Fermi-LAT excess. 

We fit to the Fermi-LAT excess energy spectrum in [3] , 
assuming the DM halo follows a generalised NFW pro¬ 
file with 7 = 1.26, r s = 20 kpc, r© = 8.5 kpc and 
Pq = 0.4 GeV cm -3 . We shower the annihilation prod¬ 
ucts with Pythia 8.186 [29]. For this model, we ob¬ 
tain todm = 44.9^4 0 GeV. As in [3], we find that 
values up to ttidm — 65 GeV provide a reasonable fit 
(p-value > 0.05). For todm = 45 GeV and the halo 
parameters mentioned, the preferred annihilation cross- 
section is ( av) = 3.2 ± 0.4 x 10~ 26 cm 3 s _1 . This is a 
factor of two larger than values in Ref. [3] because we 
assume \ is a Dirac fermion while Ref. [TTj assumed a 
Majorana fermion. 

The annihilation cross-section for \X —> A —> qq is 

/ \ _ ffDMffSM w DM _ 

1 )q 2^2 (M 2 _ 4m 2 )M )2 + M 2 aT 2 a 

( 2 ) 

This equation allows us to map (av) = J2 q ( av )q to the 
parameters in our model. The shaded blue bands in 
Fig-! show the region consistent with the Fermi-LAT 
excess. In all panels, we assumed ttidm = 45 GeV and 
(av) = 1.4 to 3.3x 10 -26 cm 3 s _1 . The lower value follows 
from variations in the halo parameters, principally pq 
which may be as large as 0.56 GeV cm -3 [33] (the anni¬ 
hilation flux $ scales as $ oc p^,M( av ))- The upper value 
follows from the Fermi-LAT 95% CL upper limit on (av) 
from dwarf spheroidal galaxies [35] . 

To compare the region consistent with the Fermi-LAT 
excess with the Mt 2 search, we establish both expected 
and observed 90% CL limits. These are given by the 
dotted black and solid red lines, respectively in Fig. [5] 
To quantify the compatibility of the expected and ob¬ 
served limits we also determine the expected ±ler and 
±2cr bands (shaded green and yellow respectively) with a 
toy experiment technique using the reported background 
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FIG. 2. Comparison of the region consistent with the Fermi-LAT excess (blue shaded) and the observed (red) and expected 
(black dotted) 90% CL exclusion contours of the Mt 2 analysis. The green and yellow shaded regions show the ±ler and ±2 ct 
bands of the expected Mt 2 limits. The black dashed and dot-dashed lines show projected limits at 13 TeV. The region above 
the lines is excluded. In each panel, we fix rriDM = 45 GeV and one of {Ma, <?dm, Ssm}, as indicated. 


uncertainties in Ref. :25j. The expected and observed 
limits also include a 20% systematic uncertainty on the 
signal yield, which is typical for these searches mug. For 
both bands we have validated our implementation with 
the Mt 2 public results and find good agreement. Based 
on the expected sensitivity of the Mt 2 search shown in 
Fig. [lj we have chosen gsM = 2 (left), (?dm = 1 (middle) 
and M a = 150 GeV (right) to illustrate the constraints 
the Mt 2 search places on the Fermi-LAT excess. A res¬ 
onance feature when Ma ~ 2?tidm is seen in both the 
Fermi-LAT region and the Mr 2 limit in the Ma — <7dm 
and Ma — <7sm planes. Outside this region, the excess is 
consistent with goM ~ gsM ~ 0(1). Owing to the off- 
shell suppression of the production cross-section, these 
searches cannot place relevant constraints on the region 
below Ma < 2todm- 

Our observed limit for the Mt 2 search is approxi¬ 
mately 2cr weaker than our expected limit. This is com¬ 
patible with [ 35] , where the observed limit for direct pro¬ 
duction of light squarks is also weaker than expected. In 
contrast, expected and observed limits are similar for the 
mono-jet analysis. This suggests that the weaker limit is 
caused by statistical fluctuations in the background es¬ 
timates in some of the phase space regions probed by 
the Mt 2 search that are inaccessible to the mono-jet 
search. Of course, we cannot exclude the possibility that 
a DM signal causes the weaker limit, but given that this 
is a 2cr effect, additional data are required to draw any 
significant conclusion. 

Even with a ~ 2cr weaker observed limit than expected, 
the Mt 2 search still excludes a significant fraction of the 
Fermi-LAT excess region for 2 todm Ma < 400 GeV. 
For gsM = 2, Mt 2 excludes all of the excess region above 
Ma = 107 GeV (left panel), while for ^dm = 1 mediator 
masses compatible with the excess above 177 GeV (mid¬ 
dle panel) are excluded. The right panel shows that Mt 2 
is able to exclude all of the excess region for g DM < 0.93 
for an illustrative mediator mass of Ma =150 GeV. In 
these panels we assumed that todm = 45 GeV but similar 


conclusions are found for values up to todm = 65 GeV. In 
fact, Fig. [T] demonstrates that the Mr 2 limits have little 
dependence on todm for todm 1$ 125 GeV. 

To illustrate how the Fermi-LAT excess parameter 
space might be probed in the future, we also provide pro¬ 
jected sensitivities of the Mt 2 search. The basis for these 
extrapolations are the 8 TeV limits, which are rescaled 
assuming that the underlying performance of the search 
in terms of signal efficiency and background suppression 
remains unchanged. These assumptions were also used 
in Ref. [16: and form the basis of Collider Reach [37]. 
Figure shows the projected limits for an early start¬ 
up scenario assuming 13 TeV and 30 ftV 1 (black dashed) 
and a long-term scenario with 13 TeV and 300 fhC 1 (black 
dot-dashed). The increase in energy and luminosity will 
enable this search to significantly increase its sensitivity. 
Assuming that search performance is maintained, it will 
be possible to probe almost all of the region Ma > 2?ttdm 
compatible with the Fermi-LAT excess. 

Finally, we see that the projected limits do not con¬ 
strain the region Ma < 2todm- This implies that this 
search will not be able to probe the ‘cascade-annihilation’ 
models that explain the Fermi-LAT excess (see e.g. [38]). 
In these models, a pair of mediating particles are pro¬ 
duced on-shell, requiring Ma < todm- 

DISCUSSION 

Although the mono-jet search is the most prominent 
search for DM at the LHC, we have shown that the 
multi-jet plus $ T search, Mt 2 > provides more strin¬ 
gent constraints on DM production for a pseudoscalar 
mediator. The additional sensitivity of the multi-jet 
search originates from binning the search into cate¬ 
gories of jet-multiplicity and kinematic variables like 
Hr and Mr 2 , as well as from extending to higher jet- 
multiplicities than the one- or two-jet final state probed 
by the mono-jet search. This is especially relevant for 
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gluon-fusion produced signal models, including the pseu¬ 
doscalar model discussed here and models with a scalar 
mediator. Powheg has the exact top-quark mass depen¬ 
dence for DM pair production and one parton in the final 
state. Topologies with higher jet-multiplicities rely on 
the parton shower for additional jets, which could intro¬ 
duce uncertainty not fully accounted by our analysis. To 
estimate the impact of the parton shower producing too 
many energetic jets, we performed a re-weighting of the 
jet-multiplicity distribution of the signal events. As an 
extreme variation, we consider the scenario where all 3- 
5 jet events are moved to the two-jet categories (while 
keeping other kinematic variables fixed) and find that 
the resulting limit is still contained in the expected ±lcr 
bands in Fig. [2] The robustness of the M T2 search against 
such variations arises from the design criteria of inclu¬ 
sive searches, which typically require that final states 
with differing jet-multiplicity have a similar sensitivity. 
Therefore, signal events assigned to the wrong category 
contribute to the analysis with a similar weight and thus 
maintain the overall performance of the search. Although 
we find that our conclusions are unchanged by this re¬ 
weighting, it would be highly desirable to have theoreti¬ 
cal tools that include the full top-quark mass dependence 
in events with multiple partons in the final state. 

The Fermi-LAT gamma-ray excess remains an enigma. 
It may be straightforwardly explained with simple DM 
models involving the exchange of an s-channel pseu¬ 
doscalar mediator, but unfortunately, the Fermi-LAT 
data are not sufficient to exclude mundane explanations 
with astrophysical sources. A DM signal in a complemen¬ 
tary experiment is required to confirm a DM origin. Con¬ 
fronting the Mt 2 limits against parameter space favoured 
by the Fermi-LAT excess shows that the LHC provides 
crucial input on pseudoscalar models. We demonstrated 
that for Ma > 2todm , much of the parameter space is al¬ 
ready constrained. Our 13 TeV projections indicate that 
essentially all of the region Ma > 2 todm will be probed 
by the next LHC run. 

As multi-jet searches consider events with at least two 
jets in the final state, the overlap with mono-jet searches, 
which allow up to two jets in the final state, is small. For 
this reason, we strongly recommend that the phase space 
covered by mono-jet and multi-jet searches is combined 
in a single search to further improve the sensitivity of the 
LHC to DM production. 
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